Entanglement Criterion for Coherent Subtraction and Coherent Addition 

Bipartite Continuous variable States 



(N 

o 

(N 

> 
O 

(N 



> 

(N 
00 

in 



X 



Li-zhen Jiang 1 , Xiao-yu Chen 1 , Tian-yu Ye 1 , Fang-yu Hong 2 

1 College of Information and Electronic Engineering, Zhejiang Gongshang University, Hangzhou, Zhejiang 310018, China 
2 Department of Physics, Zhejiang Sci-Tech University, Hangzhou, Zhejiang 310018, China 



Abstract 

Photon subtraction and addition are experimental means 
of generating non-Gaussian states from Gaussian states. 
Coherent subtraction or addition is a combination of pho- 
ton subtractions or additions. The resultant states are 
quite general non-Gaussian states. The states can be pho- 
ton number entangled states with arbitrary coefficients. 
We derive the entanglement conditions for several classes 
of coherent subtraction or coherent addition bipartite con- 
tinuous variable states. One of the entanglement condi- 
tions is necessary and sufficient. 

PACS number(s): 03.67.Mn; 03.65.Ud 



1 Introduction 

Quantum continuous variable entanglement is essential for 
quantum information tasks realized with quantum optical 
means. A great deal of contributions have been devoted 
to the entanglement criteria of Gaussian states p] [2] [3] 
[4] [5] which are continuous variable states that can be 
determined only by their first and second moments. Re- 
cently, non-Gaussian states attract much attentions for its 
good performance in quantum teleportation [B] , dense cod- 
ing [7], quantum computation [5] and nonlocality test [5]. 
Several entanglement criterion have also been proposed for 
non-Gaussian states [TU] [H] [12] [33] [HI- However, types 
of non-Gaussian states that have been studied are sporadic 
and limited by experimental conditions and theoretical in- 
terest until now. We will introduce the general ideas of co- 
herent subtraction and addition to describe non-Gaussian 
states and study their entanglement conditions. 

To generate non-Gaussian state from Gaussian state, 
one of the means is to subtract (or add) photon. The 
photon subtraction (addition) state is a™a™ pg^i ^2 
(a\ m a2 pca™a™'), where pa is a Gaussian state (here- 
after referred as Gaussian kernel) , a\ and ai are the 
annihilation operators of the two modes, a\ and a\ the 
creation operators. The replacement of a; with the 
superposition fa, + ra\ has been proposed for quan- 



tum state engineering [15] . The most generic defini- 
tion of bipartite photon coherent subtraction and ad- 
dition state is psa 



CsaPgC\ a , where C S a 



'SA' 

V°° A „1™i„tro 2 m 3 m 4 - f y. 

plex coefficients A mim2m3mi . This is so generic a defini- 
tion that all two-mode continuous variable states can be 
expressed in such form. The photon coherent subtraction 



state is p s = Csp G C\, with C s = J2 mi m 2 A mi-> 
The photon coherent addition state is pa = CaPgC\, with 



y 



m 2 ■Am 1 m 2 a l 



fmi „t™*2 



a; 



More specifically, if the 



two modes are symmetry, the symmetric coherent subtrac- 
tion and addition states are 



Ea m m II \ "> a* in in 
A m a\ a 2 p G A n a[ a 2 

m— n— 

oo oo 

Ea \m im II \ a* n « 
A m a[ a 2 p G ^A^a. 



2 > 



m=0 



n=0 



respectively. The symmetric coherent subtraction or co- 
herent addition state set is still large enough to contain 
many interesting non-Gaussian states. In this paper, we 
will study the entanglement of several symmetric coherent 
subtraction or addition states. Section 2 is to introduce 
the realignment entanglement criterion for non-Gaussian 
states, with a new form of the criterion in coherent state 
representation. Section 3 is devoted to the entanglement 
conditions of coherent subtraction state when the Gaus- 
sian kernel is a two-mode squeezed thermal state (TMST), 
when the Gaussian kernel is in its standard form. In sec- 
tion 4, we investigate the entanglement conditions of co- 
herent addition states such as photon number entangled 
states (PNES) evolved in thermal noise and amplitude 
damping environment, coherent addition of TMST. Con- 
clusions are drawn in section 5. 
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2 Realignment entanglement cri- 
terion in coherent state repre- 
sentation 

The realignment entanglement criterion has been stud- 
ied for continuous variable system recently |16j . The 
derivation is based on Fock space representation. We 
here give the criterion based on coherent state represen- 
tation. A bipartite continuous variable state p is com- 
pletely specified by its characteristic function x(p) = 
Tr[pT>(p)}, where T>{p) = exp(/i«r — p*a) is the dis- 
placement operator, and a = (ai,a 2 ) T , p = (pi,p 2 ) for 
lxl two-mode state. A Gaussian state is completely 
determined by its first and second moments. The char- 
acteristic function of a Gaussian state pa with nulli- 
fied first moments (the first moments are not relevant 
to the entanglement and can be removed by displace- 
ment) is xg(p) = ex p[ — ^(p, P*)j(p, H*) T ]i where 7 is 
the 4x4 complex covariance matrix of pa- Denote a = 
{ai,a 2 ), and denote the two mode coherent state as 
|a) = |ai,ct2)- In coherent state representation, we 
have (a\ PG \(3) = J [ilLi ^T"] Xg{p) {a\T>(-fx) \fi) = 
_i= exp [_M^I + (3) r (a*, /3) T ]. Where 

T = 01 ® I 2 + 03 ® hl'^ 1 ^ <E> h, 7' = 7 + <E> h- 
Here <n are the Pauli matrices, I 2 is the 2x2 iden- 
tity matrix. The realigned density matrix pan is the in- 
terchange of second and third subscripts of pc, namely, 
(aia 2 \ pgr I/81/J2) = {ail3\\pG \a2P2) ■ Thus 



3 Coherent subtraction 

A 1 x 1 bipartite Gaussian state can be transformed to its 
first standard form (denoted as p G ) by local operations. 
Further, it is always possible to transform p G to its second 
standard form p G by proper local squeezing operations pQ. 
To prove the necessary condition of entanglement for a 
Gaussian state, it is shown that a separable Gaussian state 
Pq takes the form of 

Pg = J P(ai,a 2 )\ai : a 2 ) (ai : a 2 \d 2 aid 2 a 2 , (1) 

where P{a\ 1 a 2 ) is a Gaussian function hence positive def- 
inite, \a>i,a 2 ) is the product of coherent states. 

Photon subtraction from Gaussian state p G will lead 
to a non-Gaussian state aiCL2p G 1 a\a 2 - More generic pho- 
ton subtraction is the coherent subtraction, or the su- 
perposition of all possible subtraction operations, namely, 
J2m m'=o A m m'aT a 2 n • The non-Gaussian state produced 
from Gaussian kernel p G is 



Ps 



E 



.4 



mm 
mm ,Uj X u 2 



y G '5>:„,4"4"'. (2) 



Then state ps is separable when p G is separable, since we 
have 



Ps 



J P(a 1 ,a 2 ) 



E 



m m 
^mm'^i ™2 



1 



(a I Pgr 



Vdet 7' 



7 exp[- 



l«| 2 + l/3| 2 



where Tr = ZTZ, with Zn = Z44 = Z 2 $ = Z$ 2 = 1 
and all the other entries of the 4x4 matrix Z are 
zeros. The realignment criterion of entanglement is 



Try PrPr > 1) from which we can derive a condition of 
entanglementTrpfl > 1. The later is easy for use but may 
detect less number of entangled states as the former. Here 
pa is the realignment of p. For a Gaussian state, we have 
the realignment entanglement sufficient condition 



Trp, 



/det 



GR 



det 7' 



> 1, 



where 7^ = er 3 ® I 2 (T R - a± ® J 2 ) V3 ® 7 2 . 

The non-Gaussian state produced by coherent subtrac- 
tion and addition from Gaussian state can be written as 
derivative of functional of Gaussian state. The method to 
derive the entanglement criterion of non-Gaussian states 
can be found in Ref.|16 . 



x \ai, a 2 ) (ax, a 2 | d 2 a\d 2 a 2 , (3) 

which is explicitly separable. Hence a non-Gaussian state 
prepared by coherent subtraction from second standard 
form Gaussian state is separable if the kernel Gaussian 
state is separable. A Gaussian state with first standard 
form p G differs from its second standard form by local 
squeezing, p G = Si{ri)S 2 {r 2 )p G S\(ri)S\(r 2 ) with some 
properly chosen squeezing parameters ri and r 2 , and Si, S 2 
are the single-mode squeezing operators. Keep in mind 
that the explicit separable decomposition (TTJ) is only for 
the Gaussian state with second standard form. So p G may 
not be diagonal in coherent state representation even when 
it is separable. 

3.1 Coherent subtraction of two mode 
squeezed thermal state 

A TMST is p st = S(r)p th S^{r), where S(r) = exp[r(a t 1 4- 
0102)] is the two mode squeezing operator, p t h = Pth,\ ® 
pth,2 is the thermal state with p t h,i — (1— Vi) J2 n v i \ n ) ( n \ ■ 
The average photon number of state pth,i is Ni — Vi/(1 — 
Vi). Consider the case of symmetric modes, we have 
N% — N 2 = N. The characteristic function of p s t is 
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Xst(p) = exp[— /i*)7(/i, /i*) T ], where the complex co- 
variance matrix 7 = bo&i <g> I 2 + C\I 2 ® <j\, with b = 
(N + ±)cosh2r, ci = -(N + ±)sinh2r. The charac- 
teristic function can also be written as x{PiiP2) == 
exp[-±(^{,/if ,^nf)M(ix[,nf , ni, fJ,^) T ], where (4 and 
p R are imaginary and real parts of pi . The real symmetric 
correlation matrix is M = 2boI 2 ® I2 — 1c\<j\ <g> 03. Notice 
that the first standard form M 1 and the second standard 
form M 11 are equal to M[T] for symmetric TMST. We 
have p\ t = p\\ = p st for symmetric TMST. The necessary 
and sufficient condition of the separability of the symmet- 
ric TMST is 

26 + 2ci>l. (4) 



Lets consider coherent subtraction state 



EA m m' A* 

m,m' Jimm U l u 2 Pst 2~in.n' I± nn' a l 



PSst — 



If {Hi is ful- 



filled, pst is separable, then ps s t is separable since it can 
be written in the form of Eq.(3). 

The inverse problem is that if ps s t is entangled or not 
when p st is entangled. We consider the case of symmet- 
ric coherent subtraction of photon, the symmetric non- 
Gaussian state is pssst = J2 m A m d^af p st J2 n A *> 



a a!, . 



The realignment entanglement criterion for unnormalized 
PSSst is [H] 

^f 0p R > OP. (5) 



where 



,0 



a 1 



3.2 Symmetric coherent subtraction of 
standard form Gaussian state 

The first standard form Gaussian state p G is com- 
pletely characterized by its complex covariance matrix 
7 = b <Ji ® I 2 + C\I 2 <S> <J\ + c 2 <7\ ® (Ji, where b = 
Tr[p G {a\ai + |)] = Tr[p G (a\a 2 + 5)], c x = Tr\p G a x a 2 \, 
c 2 = Tr[pQa\a 2 ]. The real symmetric correlation matrix 
is M 1 — 2b I 2 <S> J2 — 2cicri ® 0-3 + 2c 2 ai ® ^2- The neces- 
sary and sufficient condition of the entanglement of p G is 



t = 777 1 — t > 1. For symmetric coherent subtraction 

state pss = XL^m a r ia 2™PGL„ A* a i™4"' we consider if 
it is separable at the r = 1 and C2 ^ 0, the boundary 
of separability of its Gaussian kernel. The entanglement 
criterion is 

OP R > OP, (6) 

where P R = exp[b(r) 1 rj 2 +£i£ y2 )~ci(r] 1 £ 1 +ri 2 £ 2 ) + c 2 (r] 1 £ 2 + 
rfrZi)], P = exp[-ci(77iry 2 +Ci6)+ & ( ? yiCi+ ? 726)+C2(f7i??2- 
£162)]) with b = b a 



§. Then 



i(m,n) 



0P K = Y^m\n\A m A* n J2 



fe=0 



k l° 2 



min(m,n)— fc 
i=0 



m—k\ I 11 — k \ , n +m-2k-2l 21 



Emin(m' ,n') 
2=0 



( a .2j_ a .m'+n'-2i)_ Notice that 



4(6o+ci) 2 

, P R = eML(viti + 

f?2&) + ^(?7i»72 + P = cxp[(fe - §)(»7i£i + %6) - 

ci(r?ir? 2 + with if = ±[(& + c x )r + (6 - c x ) - 1], 

L = 2 [(bo — c\) — (bo + ci)t). Then The entanglement con- 
dition is 26o + 2ci < 1, so we keep 60 — c\ invariant. Denote 
e = 1 — 26o — 2co, then e is a small positive quantity for 
entangled symmetric TMST near the separable boundary. 
Expanding all the quantities in ([?]) to the first order of e, 
we have bo + co = \ — e, r«l + 4e, K = z+ L = z— |, 
bo — I = z — |,ci = — z — ~, where z = |(6q — ci — |). 
Notice that up to the first order of e, we have OP R = OP, 
while \j~r > 1. Hence (O is fulfilled even for infinitive 
small e. Thus pssst is entangled if p s t is entangled. We 
have proven that for any symmetric coherent subtraction 
state pssst = E m A m a?a 2 n p st En A >\ n 4 n produced where 
from symmetric TMST, the necessary and sufficient of 
separability is (01 . 



The same expression can be found for OP with the in 
terchange of b and —c\. Consider the function f m 'n'(x) = 
m! \ I n 

1 ^ l 

fm'm'(x) = 0, we only need to consider ml 7^ n'. For odd 
ml + n' , one of the roots for equation f m 'n'ip) = is 
1=1; For even ml + n', x = ±1 are the roots of equation 
fm'n'(x) — 0. We have verified numerically that there are 
no other real roots for all m',n' < 200. Since /,„'„< (0) = 1, 
we have f m > n '(x) > for all x € (—1,1). From r = 1, 

we have b + c\ - 



\ > 0, thus -f G (0,1) as 



ci is assumed to be negative (for positive c\, the realign- 
ment criterion Trpn > 1 should be modified [H]), then 
fm'n'(— %) > 0. At r = 1, the realignment entanglement 
criterion OP fl — OP > can be rewritten as 



min(m.n) 
m>n fe=0 



?7J 



^. „2kin+m-2k f ( C l \ ^ r, 

XC 2 Jm-k,n-k(--^) > U, 



= m!n!(j4 m A* + AnA^). 



(7) 



(8) 



If B„ m > for all m, n and -B mn > at least for one 
pair of m,n (m ^ n), the coherent subtracted state pss. 
One of the case is that the real and/or imaginary parts 
of the coefficients A m are non negative and at least two 
of A m have positive real and/or imaginary parts. Coher- 
ent subtraction means at least two of the coefficients A m 
are non zeros. We conclude that if the coefficients have 
positive real and/or imaginary parts, the symmetric non- 
Gaussian state produced from the coherent subtraction of 
an edge separable first standard form of Gaussian state is 
entangled. 
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4 Coherent addition 

Unlike the coherent subtraction state, the sufficient con- 
dition of separability for a coherent addition state is dif- 
ficult to be found even when the Gaussian kernel pc is 
in its second standard form. We will consider the suf- 
ficient condition of entanglement for symmetric coherent 

tm„lTO n „V°° A*n n n n 
2 PG l_in=Q A n a \ a 2 



addition state psa = Em=o ^™ a i a 
in the following. An initially prepared non-Gaussian en- 
tangled state will become separable after it interact with 
the thermal noise and amplitude damping environment. 
The evolution of the state could be specified by the evo- 
lution of its characteristic function [18]. The realignment 
entanglement criterion for time evolution state is [16] 



det(V « ) -Oexp[-ivfiv T ] > Oexp[-iv 7 'v T ] 



det( 7t ') 



(9) 



where v =(e x , s 2 , -Ci, ~Ca) 

^ ID I 



o 



<i=0 



= E 

n = 7 



, A m A* n O mn yith 



i t - Y Z'i m Z'i t -^)i with i t = e -W + (n+l)(l 



I 2 , Y Rt = {Z'i^Z 1 + I 2 ® a 1 + o-i <8> /2)- 1 and Z' - 
(0-3 (g> 72)-^(c3 ® 12)- Here n is the thermal noise, the am- 
plitude damping. coefficient is merged in the unit of time 
t. 

4.1 Photon number entangled state 

A special case of symmetric coherent addition state is 
the PNES when the Gaussian kernel pa is the vac- 
uum state. We have ppnes = \ip) (V 1 ! where \ip) = 



m mm) 



For vacuum Gaussian kernel, 



7 = 



0-1(8/2, 7t = n t oi®I 2 , i m = z?^[{l-n t )h®ai+nt<Jx® 
/ 2 ], where ra t = n(l - e _t ) + 1. Then f2 = -zl 2 ® <T\ + (1- 
z)<7i ®/ 2 , with 2 = e- t /{2n t - 1) = e _ 7[2n(l - e"') + 1]. 
Denote F = -±vOv T = z( £ i£ 2 + C1C2) + (1 - + 
£2C2), then 



(9, 



1 



(m + n)\ 
mini 



min(m,n) 
k=0 




2 4 
N 



0.2 0.4 0.6 0.8 1.0 

IcJ 2 



Figure 1: (a)The critical damping time of separability for 
TMC, N is the average energy of input state, it is a little 
bit less than A. From up to bottom, the thermal noise 
n=10< - 5), 10< - 4), 10< - 3), 10< - 2), 10< - 1), respectively. 
(b)The critical damping time of separability for photonic 
Bell-state. The solid line is the realignment criterion, the 
dash line is the Simon's criterion. 



For a two mode squeezed thermal state evolving in 
amplitude damping and thermal noise channel, the ini- 
tial state is \ip) = y/l — A 2 Em=o ^ m \ mm ) 7 tnus C m = 
vl — A 2 A" 1 . By interchanging the order of summations in 
([TO]) , the entanglement criterion of the final state reduces 
to 

(id 



(2n t -l)(l + A-2Az) 



> 1. 



Condition (lll[) can also be derived from the well know en- 
tanglement condition of Gaussian state, for the final state 
is Gaussian. This example shows that the choice of Gaus- 
sian kernel may not be relevant for the entanglement cri- 
terion. 

For a two mode coherently correlated state (TMC) with 
Poisson coefficients evolves in thermal noise and ampli- 
tude damping channel, the unnormalized initial state is 



W = Em=o T^T \mm) , then C m = 



s(A)m! 



where g(x) 



(1 - zf k z 



2k n+m-2k 



We also get O mn exp[-±v 7 'v T ] = m).n\5 mn , y ^"*) = 
„ 1 , . The realignment entanglement criterion for time 

2n t — 1 & ° 

evolution PNES is 



Em(^rf) 2 - The function g(x) converges fast than expo- 
nential function. The entanglement criterion (fit)]) leads 
to 

e 2Xz g(X(l-z))>(2n t -l)g(X). (12) 

The critical separable time is shown in Fig. 1(a) for differ- 
ent noise n. Comparing with the direct numerical calcula- 
■tion [17] , we confirm that the numerical calculation of the 
realignment entanglement criterion of TMC is valid. 

Another example is the evolution of photonic Bell-state 
\^b) = Co 1 00) + ci 1 11) in thermal noise and amplitude 
damping environment. The entanglement criterion (|10p 
leads to 



r , min(ro,n) 



2m - 1 ^ 



k=0 



where C m = 



2\A„ 



(l-z) 



2k n-\-m—2k 



> 1. 



(10) 



I c 1 + 2Re(c cl)z + |ci| 2 (1 - 2z + 2z 2 ) > 2n t - 1. (13) 

The critical separable time is shown in Fig. 1(b), where the 
Simon's criterion is also shown for comparison [2]. The 
similar curve has been obtained by Nha's criterion [14]. 
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4.2 Coherent addition of 
squeezed thermal state 



two mode 



The coherent addition of symmetric TMST is 
PASst = Em=o A m a\ m al m p st £^° =0 A* a? c$ . The re- 
alignment entanglement criterion is still ([9]), where 



A A* 



■i[(6o 

/ 2 + 



= =o O = i[(6 -ci) 

+ ci)r + (& - ci) 
C1/2 <8> ci, with t 



0"1, 



(6 + ci)r]t7i ® 7 2 

7' = (&0 + j) ! & x 2 -T cij 2 W uii wiuii 7 = 4(b 0+Cl )2 

We consider the case that the kernel is a symmet- 
ric TMST of critical separable, that is t = 1. Then 



= — c\&\ ® h ~ (bo 
criterion © then reads 



cr\. The entanglement 



E 



Bmn (bo 



1 



n+m 



/"m , n ( ~ 



Cl 



&0 



>0, 



1. 



where i? m „ is defined as in ©. Notice that when r 
we have ~ b +1 = -^+1 e (Oj 1)> nere c i < is assumed 
as before. Hence f m ,n(— b ^_i ) > as we have argued. So 
if all B mn are nonnegative (at least one of them is nonzero 
thus it is positive), the state should be entangled. One of 
the cases is that if only two of the coefficients are nonzero, 
say A m and A n , the entanglement condition will be 

A m A n + A n A m > 0. 



Denote A m = \A m \e 
duces to 

cos(</9 



, the entanglement condition re- 
<p n ) > 0. 



for any m 7^ n. Thus if the argument difference of the 
two complex coefficients is an acute angle, the state is 
entangled. Borrowing term from optics, we can say that 
if the coherent addition constructively interferes, the state 
is entangled. 

5 Conclusions 

We have introduced the ideas of coherent subtraction and 
coherent addition to describe many kinds of bipartite non- 
Gaussian states and investigate their inseparability. The 
non-Gaussian states produced from their Gaussian kernels 
are quite generic. The results of the paper are of several 
folds. For any coherent subtraction state produced form 
the second standard form [1] Gaussian state, we prove 
that it is always separable when the Gaussian kernel is 
separable. For any symmetric coherent subtraction state 
produced from the symmetric two-mode squeezed thermal 
state, the necessary and sufficient separable condition is 
obtained and it is exactly the same as that of its kernel 
state. A symmetric coherent subtraction state produced 
from the symmetric first standard form Gaussian state (ex- 
cept two-mode squeezed thermal state) is entangled when 
the Gaussian kernel is at the boundary of separability and 



all the superposition coefficients have positive real and/or 
imaginary parts. Thus a proper arranged coherent sub- 
traction produces entanglement from a separable state. 
For photon number entangled state evolving in thermal 
noise and amplitude damping channel, we give the analyt- 
ical entanglement condition. As applications, we present 
the separable conditions of two mode coherently correlated 
states with Poisson coefficients and photonic Bell states 
evolving in thermal noise and amplitude damping envi- 
ronment. Any two term coherent addition state produced 
form boundary separable symmetric two-mode squeezed 
thermal state is entangled if the two complex coefficients 
give rise to constructive interference. Further works are 
desirable for other non-Gaussian states or using entangle- 
ment criterion other than realignment criterion. 
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